function path=ouprocess(N)

m=1; alpha=1; y0=2; beta=1;
time=linspace(0,1,N);
deltat=1/(N-1);
path=zeros(size(time));
path(1)=y0;
hold on
for pathnum=1:100
for ii=2:N
    path(ii)=path(ii-1)+alpha*(m-path(ii-1))*deltat+beta*randn(1)*sqrt(deltat);
end
plot(time,path);
end
